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In recent experiments, localized and stationary pulses have been generated in second-order nonlin- 
ear processes with femtosecond pulses, whose asymptotic features relate with those of nondiffracting 
and nondispersing polychromatic Bessel beams in linear dispersive media. We investigate on the 
nature of these linear waves, and show that they can be identified with the X-shaped (0-shaped) 
modes of the hyperbolic (elliptic) wave equation in media with normal (anomalous) dispersion. De- 
pending on the relative strengths of mode phase mismatch, group velocity mismatch with respect 
to a plane pulse, and of the defeated group velocity dispersion, these modes can adopt the form of 
pulsed Bessel beams, focus wave modes, and X- waves (0-waves), respectively. 

PACS numbers: 42.65.Re, 42.65.Tg 



I. INTRODUCTION 



Stationary, temporally and spatially localized, X- 
shaped optical wave packets, having a duration of a few 
tens of femtoseconds and spot size of a few microns, have 
been recently observed to be spontaneously generated in 
dispersive nonlinear materials from a standard laser wave 
packet 0,1111. Balancing between second-order or Kerr 
nonlinearity, group velocity dispersion (GVD), angular 
dispersion and diffraction, has been suggested to act as 
a kind of mode-locking mechanism that drives pulse re- 
shaping and keeps the interacting waves trapped, phase 
and group-matched 0, |5i,i6s|. 

The purpose of the present paper is to investigate on 
the nature of these waves. The main hypothesis underly- 
ing our investigation is that these nonlinearly generated 
X-shaped waves behave asymptotically as linear waves. 
This assumption is based, first, on the observed station- 
arity, not only of the central hump of the wave packet, 
but also of its asymptotic, low-intensity, conical part 
stationarity that cannot be attributed to nonlin- 
ear wave interactions, but to some linear mechanism of 
compensation between material and angular dispersion. 
Indeed, several kinds of linear polychromatic versions of 
Bessel bea ms B , as Bessel-X pulses 8, 9], pulsed Bessel 
beams Rfl [ll|, subcycle Bessel-X pulses or focus wave 
modes |Ti|, and envelope X waves [13, with the capa- 
bility of maintaining transversal and temporal (longitu- 
dinal) localization in dispersive linear media, have been 
described in recent years (for an unified description and 
the extension to media with anomalous dispersion, see 
also Ref. 0). In contrast to free-space X- waves and 
Bessel-X pulses 0, stationarity in dispersive me- 
dia requires the introduction of an appropriate amount 
of cone angle dispersion that leads to the cancellation of 
material GVD with cone angle dispersion-induced GVD 
H H [13, m d El m. Second, polychromatic Bessel 
beams, with or without angular dispersion E3|i have the 



ability of propagating at rather arbitrary effective phase 
and group velocities in dispersive media, as has to be 
done by the phase matched and mutually trapped funda- 
mental 3] and second harmonic [^ nonlinearly generated 
X waves. 

For these reasons, in this paper we present a new and 
more comprehensive description of localized and station- 
ary optical waves in linear dispersive media, henceforth 
called wave modes, that is particularly suitable for un- 
derstanding and predicting the spatiotemporal features 
of the nonlinear X waves generated in experiments. On 
the linear hand, this description allows us to predict the 
existence of new kinds of wave modes, and classify all 
them according to the values of a few physically mean- 
ingful parameters. 

Each wave mode is specified by the values of the de- 
feated material GVD, the mode group velocity mismatch 
(GVM) and phase mismatch (PM) with respect to a 
plane pulse of the same carrier frequency in the same 
medium. Wave modes are then shown (Section to 
belong to two broad categories: hyperbolic modes, with 
X-shaped spatiotemporal structure, if material dispersion 
is normal, or elliptic modes, with 0-shaped structure, if 
material dispersion is anomalous 18]. In Section UTTl we 
show that each wave mode can adopt the approximate 
form of 1) a pulsed Bessel beam (PBB), 2) an envelope 
focus wave mode (ePWM), or 3) an envelope X (eX) wave 
in normally dispersive media [envelope O (eO) wave in 
anomalously dispersive media] , according that the mode 
bandwidth makes PM, GVM or defeated GVD, respec- 
tively, to be dominant mode characteristic on propaga- 
tion. This classification allows us to understand the spa- 
tiotemporal features of wave modes in dispersive media in 
terms of a few parameters (the characteristic PM, GVM 
and GVD lengths), including modes with mixed pulsed 
Bessel, focus wave mode, and X-like (0-like) structure. 

The above description is obtained from the paraxial 
approximation to wave propagation. We choose this ap- 
proach because of its wider use in nonlinear optics, and 
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because it leads to simpler expressions in terms of param- 
eters directly linked to the physically relevant properties 
of the mode and dispersive medium. In Section IIVI we 
compare the paraxial and the more exact nonparaxial 
approaches, to show that the paraxial approach is accu- 
rate enough for the description of wave modes currently 
generated by linear optical devices and in nonlin- 

ear wave mixing processes 1, 2, 3]. 

II. WAVE MODES OF THE PARAXIAL WAVE 
EQUATION 

We start by considering the propagation of 
a three-dimensional wave packet E{x±,z,t) — 
A{x±, z,t) eTq){—iuJot + ikoz) [x± = {x,y)] of a certain 
optical carrier frequency loq, subject to the effects of 
diffraction and dispersion of the material medium. 
Within the paraxial approximation, and up to second 
order in dispersion, the propagation of narrow-band 
pulses is ruled by the equation 

^^A^^^A^A-^'f^lA, (1) 

where z is the propagation direction, t — t — k'^z is the 
local time, = + dy, and k^^ = 9^ 'fc(ijj)|j^„ , with 
k{uj) the propagation constant in the medium. Eq. |^ is 
valid for a narrow envelope spectrum A{x^^ z, 51) around 
51 = w — Wo = 0, that is, for bandwidths 

All < Wo, (2) 

a condition that requires at least few carrier oscillations 
to fall within the envelope A. 

We search for stationary and localized solutions of Eq. 
([T|) in the wide sense that the intensity does not depend 
on z in a reference frame moving at some velocity. These 
solutions must then be of the form 

A[x, y, r, z) — $(a;, y,T + az) exp(— i/3z). (3) 

The free parameters a and P are assumed to be small in 
the sense that 

|a| « ki (4) 
\P\ « ko, (5) 

so that the group velocity l/(fcQ — a) and phase velocity 
i^o/{ko — 0) of the wave differ slightly from those of a 
plane pulse of the same carrier frequency in the same 
material, l/fep and ujo/ko, respectively. 

Under the assumption of asymptotic linear behavior of 
nonlinear X- waves, we can get some insight on the pos- 
sible values of a and /? of nonlinear X- waves on the only 
basis of the linear dispersive properties of the medium. 
If, for instance, a pulse of frequency up generates a sta- 
tionary and localized second harmonic pulse (wq = 2ujf) 
travelling at the same group and phase velocities as the 



fundamental pulse Q, we must have k'^ — a = k'p and 
ki)—(3 = 2kF, that is, a = kp — k'^ and /3 — Afc = kQ~2kp. 
For illustration. Fig. ^ shows the values of a and (3 of 
the second harmonic pulse in lithium triborate (LBO) as 
a function of its carrier frequency uq. Note also that \a\ 
and 1/3 1 satisfy the conditions (@J| and (O for any carrier 
frequency in the entire visible range and beyond. 
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FIG. 1: Values of a and /3 of the localized and stationary, 
second-harmonic waves of different carrier frequencies ujo for 
phase and group matching with the fundamental wave in the 
process of oo-e second harmonic generation in LBO at room 
temperature. Dispersion formulas for the refraction index are 
taken from Ref. pi)| 

In Section [l VI a nonparaxial approach to the problem 
stated above will be performed. It will be shown that 
the paraxial and nonparaxial descriptions yield substan- 
tially the same results if conditions (|2Jl, (01 and ||SJ) are 
satisfied, as is the case of the experiments and numerical 
simulations demonstrating the spontaneous generation of 
X-type waves jlSli. 

Equation with ansatz (jSJ yields 

Aj_^ - /co^o^r* + 2ikoadr<P + 2kof3^ = 0, (6) 

for the reduced envelope or the Helmholtz-type equa- 
tion A±<t + i^^(Sl)$ — 0, for its temporal spectrum 
<i?(x, y, fl), where 



K{fl) = ^ 2ko \^P + an + -k'f^n^j (7) 

will be referred to as the (transversal) dispersion relation 
since it relates the modulus K of the transversal com- 
ponent of the wave vector with the detuning 51 of each 
monochromatic wave component from the carrier fre- 
quency ajQ. For 51 such that if (51) is real, the Helmholtz 
equation admits the bounded, cylindrically symmetric, 
Bessel-type solution <I>(r, 51) — f{^l)jQ[K{^l)r], where 
/(51) is an arbitrary spectral amplitude and Jo(') the 
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Bessel function of zero order and first class |20j. By in- 
verse Fourier transform we can write the expression 



$a,/3(r,T + az) 



1 



271" JK{n) real 
X Jo[K{n)r]exp[~in{T + az)] (8) 



for the reduced envelope of the cylindrically symmetric 
wave modes, or localized, propagation invariant solutions 
of the paraxial wave equation, in the sense explained 
above. As indicated, the integration domain extends over 
frequencies such that the dispersion curve K{^1) is real. 
According to Eq. (jSjl, a wave mode ^a.js is composed of 
locked monochromatic Bessel beams whose frequencies 
and radial wave vectors are linked by a specific dispersion 
relation K(Q), and whose relative weights are determined 
by a certain spectral amplitude /(Jl). 



normal dispersion 
( > ) 




anomalous dispersion 
(V<0) 




{kg > 0), K{n) is in fact a single-branch vertical hyper- 
bola if /3 > Q!^/2fcQ, and a two-branch horizontal hyper- 
bola if /? < o? jlkg [see Fig. EJa)]. For anomalous disper- 
sion (fep < 0), K(p^ takes real values only if /? > /2k'^, 
in which case the dispersion curve is an ellipse [see Fig. 
^2p\ ■ It is also convenient to introduce the (real or imagi- 
nary) frequency gap 
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and radial wavevector gap 



Kg — J —k^k'^Vl'^. 



(9) 



(10) 



When Vlg and Kg are real, they represent actual fre- 
quency and radial wavevector gaps in the dispersion 
curve K(Vl), as illustrated in Fig. |21 In any case, their 
moduli characterize the scales of variation of the fre- 
quency and radial wavevector in the dispersion curves. 

Closely connected with the dispersion curve are the so- 
called impulse response wave modes ^')^fj{i',T + olz), or 
modes with f{Vl) — \. As seen in Fig. [3 the structure 

of in space and time closely resembles that of the 
dispersion curve in the K-VL plane, but at radial and 
temporal scales of variation determined by the reciprocal 
quantities |i^g|~^ and |^lg|~^, respectively. Eq. © with 
/(SI) = 1 and the change Vt' = ri + a/fep yields 



1 



dVL' 



27r JK{n') real 
X Jo [K{n')r] exp[-ir2'(T + az)] 
a 



X exp 



^0 



(11) 



where 



FIG. 2: Dispersion curve of the wave modes in a medium 
with (a) mormal dispersion, (b) anomalous dispersion. 



K{n') = Jkk'^ 



(12) 



As shown in Figs. Wia) and (b), the form of the dis- 
persion curve K{^) reflects the underlying hyperbolic or 
elliptic geometries of the paraxial wave equation in 
the respective cases of propagation in media with nor- 
mal or anomalous dispersion. For normal dispersion 



The integral in Eq. can be performed in all pos- 

sible cases from formulae 6.677.3 (for k'^ > 0, /3 > 
Q!V2fc(;),6.677.2 (for k'^ > 0, /3 < /2k'^) and 6.677.6 
(for k'^ < 0, /3 > a^l2k'^) of Ref. to yield the closed 
form expression for impulse response modes 



^%,{r,T + az) = — 
X exp 



1 



y/kok/^^r^ ~{t + az)' 



la , , 
— {T + az) 



^expji sj^-^sl kok'^r^ - {t + azf | 



^ C. C. 



(13) 



or, in terms of the frequency and radial wavevector gaps 

_ 1 



Q,/3 



27r 



exp(^^_^^_ C. 
iR 



exp 



la , , 
— (j^az) 



(14) 



where R = [K^r^ + flUr + azf^l'^. 
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FIG. 3; Gray-scale plot of the amplitude I'l'^^l of the impulse response wave modes, (a) Normal dispersion k'o > with 
imaginary, transversal wave vector gap Kg real), (b) Normal dispersion fco > with fi < /2k'(^ (detuning 
gap Q,g real, Kg imaginary), (c) Anomalous dispersion /cq < with (3 > /2k'i^ {Q,g and Kg real). Normalized local time and 
radial coordinate are defined as a = I^^sKt + oiz) and p = l-ft^gk, respectively 



As shown in Fig. EJa), for fc^' > and /3 > /2k'^ {Vlg 
imaginary and Kg real) , the impulse response wave mode 
is singular in the cone r ~ |(t + az)\/ ^Jkok^, is zero for 
r < |(t + az)\/ ^Jk^k'f^ (within the cone), and decays as 

1/r for r > \{t + az)\/ -^/fco^o ('^^^ '^^ ^^'^ cone). The 
radial beatings in this region, of period 2TT/Kg, are a 
consequence of the radial wave vector gap Kg . 

Figure|2Ib) shows the impulse response mode for /cq > 
and P < a^/2fcg {fig real and Kg imaginary). As in the 
previous case, the mode is singular at the cone r — | (r + 
<^z)\/\/kokQ, but damped oscillations are now temporal, 
of period 27r/f2g, as corresponds to the frequency gap 
in the dispersion curve. Out of the cone [r > |(t + 
az)\/ y/kak'^] the mode is exponentially localized. 

Modes in media with anomalous dispersion, i.e., with 
kQ < and /3 > (real ^g and Kg), exhibit rather 

different characteristics [Fig. OJc)]. These modes are no 
longer singular and of X-type, but regular and, say, of 
0-type. The damped oscillations decay temporally and 
radially as 1/t and 1/r, respectively, with periods 27r/r2g 
and 2tt / Kg. The absence of singularities is a consequence 
of the actual limitation that the elliptic dispersion curve 
imposes to the uniform spectrum /(Jl) = 1. 



III. CLASSIFICATION OF WAVE MODES 

Numerical integration of Eq. |(SJ) with a given disper- 
sion curve (specified by the values of a, (3 and fcg ) but 
different (bell-shaped) spectral amplitude functions /(fi) 
having also different (but finite) bandwidths Af7 [alter- 
natively, numerical integration of 

/•OO 

*a,/3(?-,'r + az) = / dCT$^'y r, T + az - cr) (15) 



where /(r) is the inverse Fourier transform of /(r2)], 
shows much richer and complex spatiotemporal features 
in comparison with the case of infinite bandwidth. These 
features strongly depend on the choice of the spectral 
bandwidth Af2, while no essentially new properties arise 
from the specific choice of /(il) (Gaussian, Lorentzial, 
two-side exponential. . . ). Modes with finite bandwidth 
may exhibit mixed, more or less pronounced radial and 
temporal oscillations, along with incipient or strong X- 
wave (0-wave), focus wave mode or Bessel structure, as 
explained throughout this section (see also the follow- 
ing figures). The purpose of this section is to perform 
a simple, comprehensive classification of wave modes in 
dispersive media. In the remainder of this paper, Afi 
will refer to any suitable definition of half-width of the 
bell-shaped spectral amplitude function f{^). 

Given a mode of parameters a and (3 satisfying con- 
ditions Q and ISJ, propagating in a dispersive material 
with GVD fcg , and some spectral bandwidth ASl satisfy- 
ing we have found it convenient to define the three 
following characteristic lengths: 1) the mode PM length 



Lr, 



1 



2) the mode walk-off, or GVM length 



(16) 



(17) 



measuring, respectively, the axial distances at which the 
mode becomes phase mismatched and walks off with re- 
spect to a plane pulse of the same spectrum in the same 
medium, and 3) the GVD length 



1 



k'^{Any^ 



(18) 



or distance at which the mode (invariable) duration dif- 
fers significantly from that of the (broadening) plane 
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FIG. 4: Characteristic lengths of extraordinary second- 
harmonic wave modes of different frequencies in the visible 
range that travel at the phase and group velocities of the 
ordinary fundamental waves in LBO at room temperature. 
Mode bandwidths are AQ,/ujo = 1/2ttN, with (a) A'' = 10 
and (b) = 1. 

pulse. Note that, as defined, Lp, and Ld can be 
positive or negative. In terms of the mode lengths the 
transversal dispersion relation ^ takes the form 



K{n) = y 2fco yLp' + L^'^n + j ' (19) 

where r2„ ~ fi/Ari is the normalized detuning, which 
ranges in [—1, +1] for Vl within the bandwidth Ail. Then, 
they are the values of the mode lengths Lp, and 
that determine the form of the dispersion curve within 
the spectral bandwidth, and hence the parameters that 
determine the spatiotemporal structure of the mode, as 
shown throughout this section. We analyze here three 
extreme cases, namely, 

\Lp\ <^ \Li^\,\Ld\ PM-dominated case 
\Lw\ <^ \Lp\,\Ld\ GVM-dominated case 
\Ld\ \Lp\,\Liu\ GVD-dominated case 

that represent three well-defined, opposite experimental 
situations, and that allow us also to understand, at least 
qualitatively, the features of general, intermediate cases. 

For illustration, we have evaluated the characteristic 
lengths of wave modes of different frequencies luq that 
propagate in LBO at the phase and group velocities of 
the corresponding fundamental waves of half-frequency. 
In Figs. 01 the bandwidths Ail — ojo/iirN corre- 
spond to "TV-cycle" pulses [duration ~ (Af])"^ = NTq, 
Tq — 27r/wo period] at each frequency uiq. The value 

= 10 in Fig. ^a) leads to a pulse duration (Afl)~^ ~ 
20 fs at ujQ = 3.55 fs~"^ (A = 0.53 /.tm), of the same 
order as in previous experiments and numerical simu- 
lations. Fig. 01b) shows, in contrast, the extreme case of 
"single-cycle" wave modes. Generally speaking, modes 



of long enough duration belong to, or participate mostly 
of, the PM-dominated case [as in Fig0fa) for most fre- 
quencies], modes of some (still unspecified) intermediate 
duration belong to the GVM-dominated case, and ex- 
tremely short modes to the GVD-dominated case, since 
Lp is independent on bandwidth, but Lp and Ld are in- 
versely proportional to Ail and Ail^, respectively. De- 
pending, however, on the relative values of a, /3 and /cq 
(particularly when one or two of them are very small), the 
GVM-dominated case, even the PM-dominated case, can 
extend down to the single-cycle regime [as in Fig0Jb) 
for most frequencies], or, on the contrary, the GVM- 
dominated case, even the GVD-dominated case, apply 
to considerably long modes [as in the vicinity of the two 
singularities of the Lp-curve of Fig. EJa)] . 



A. Phase-mismatch-dominated case: Pulsed Bessel 
beam type modes 

Consider first modes with |Lp| <C When 
Lp > 0, the dispersion curve within the spectral band- 
width can be approached by the real constant value 
K{n) ^ [2k^Lp^f/\ or, 

K{n)^^2k^ (if/3>0), (20) 

[see Fig. El^a)] regardless the exact dispersion curve is 
an actual hyperbola or ellipse [as in Fig. Elb)], that is, 
independently of the sign of material group velocity dis- 
persion. Wave modes under these conditions can only 
have superluminal phase velocity (/3 > 0), but super- or 
subluminal group velocity (a > or a < 0, respectively), 
and will adopt, from Eqs. ||SJ) and H20|l . the approximate 
factorized form 

<^c,AT,T + az)'^ f{T + az)Jo[y2k^r^ (21) 

of a PBB of transversal size of the order of (2fco/3)^^^^. 

Figure [SJc) shows the prototype PBB of this kind of 
wave modes [Eq. (I21ll ] with a Gaussian spectrum /(fl), 
that is, the limiting case \Lp/L.uj\ — 0, \Lp/Ld\ = 0, 
or horizontal thick lines of Figs. E^a) and (b). In 
Fig. |5fd) we show, for comparison, the wave mode with 
\Lp/Lw\ = 0.25, \Lp/Ld\ = 0.25 and with the same Gaus- 
sian spectrum, obtained numerically from Eq. ||SJ). We 
see that the wave mode preserves a spatiotemporal struc- 
ture similar to that of the prototype PBB of Fig. I3c), 
even if |Lp| is not much smaller, but simply smaller than 
|L^| and \Ld\. Small differences can be understood as 
incipient focus wave mode and 0-wave type behavior, as 
described in the following sections. 



B. Group-velocity-mismatch-dominated case: 
Envelope focus wave modes 
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FIG. 5: (a) Dispersion curve within the bandwidth for |_Lp|/|Lu,| 0, \Lp\l\Ld\ (thick curve), and for Lp/L^u = —0.25, 
Lp/Ld ~ —0.25 (thin curve), (b) The same as in (a) but also outside the bandwidth of the Gaussian spectrum (in arbitrary 
units) /(SI) = exp[— (il/An)^]. (c) and (d) Gray-scale plots of the amplitude I'&a./sl of (c) the PBB of Eq. I)21|l with spectrum 
f{n) = exp[-(n/Ar2)^] (i.e., /(r) oc exp[-(2Anr)^]) and (d) of the mode with Lp/L^ = -0.25, Lp/La = -0.25 and same 
spectrum as in (c), numerically calculated from Eq. (|5J. Normalized coordinates are a — {t + az)AQ,, p — r/ro, with 
ro = (2fco/3)-'/'. 




FIG. 6: (a) Dispersion curve within the bandwidth for Lw = 10/fco, L-w/Lp — » 0, L^/Ld (thick curve), for Lw ~ 10/feo, 
Lm/Lp = 1/8, L^/Ld ~ 1/8 (thin curve, label 1), and for = lOfco, L^/Lp = 1/3, L^jLd ~ 1/3 (thin curve, label 2). (b) 
The same as in (a) but also outside the bandwidth of the Gaussian spectrum (in arbitrary units) f{fl) — exp[— (f2/Af2)^]. 
(c) and (d) Gray-scale plots of the amplitude |3>c«,/3| of (c) the prototype eFWM [thick dispersion curve in (a)] with spectrum 
/(fl) = exp[—{fl/AQ)^], and of (d) of the mode with L^/Lp — 1/3, L^/Ld = 1/3 [thin dispersion curve 2 in (a)], numerically 
calculated from Eq. ||5J. Normalized coordinates are a = {t + Qz)Af2, p = r/ro, with ro = (2//coAf2|Ql)^/'^. 



The case ^ |Lp|, |Lc(| leads to a new kind of wave 
modes that has not been reported. The dispersion curve 
within the bandwidth is now of the form of the horizontal 
parabola K{n) ~ (2/coiu;^f^n)^^^ with vertex at = 0, 



or, 

K{n) ~ ^j2koaVL, (22) 
[see Fig. EJa)], regardless material dispersion is normal 
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[as in Fig. Efb)] or anomalous. For modes with superlu- 
minal group velocity (a > 0), the horizontal parabola 
is right-handed [as in Figs. Ela) and (b)], and left- 
handed for subluminal modes (a < 0). Independently 
of the group velocity, phase velocity can be superlumi- 
nal (/3 > 0) or subluminal {(3 < 0). In any case, their 
spatiotemporal form can be approached by Eq. ijSJ) with 
K{^l) given by Eq. (|22|) . Moreover, with the two-sided 
exponential spectrum f{fl) ~ {2tt/ Aft) cxp{—\n\/ Aft), 
Eq. © yields 



^a,f3ir,T + az) 



-ZTq 



az 



■ exp 



ikQ\a\r'^ 



2(r- 



iTo) 



(23) 

for superluminal modes (a > 0), and the complex conju- 
gate of the r.h.s. of Eq. H23|l for subluminal modes {a < 
0). In Eq. I|23|) . tq = (Aft)^^ characterizes the mode du- 
ration. The mode spot size at pulse center (t + az ^ 0) 
can be characterized by rp = (2/fcoAr2|Q;|)^/^. 

The functional form of the reduced envelope in Eq. 
(|23|l is similar to the fundamental Brittigham-Ziolkowski 
focus wave mode (FWM) 0, and as such will be 
called envelope focus wave mode (eFWM). There are, 
however, important physical differences between them, 
which can be understood for the respective expressions 
of the complete fields E of both kind of waves, namely. 



-ITq 



■ exp 



ikolalr'^ 



2(7 



X exp{—il3z) eiq){—iujot + ikoz), (24) 



for the envelope focus wave mode, 



E{r,z,t) 



-ITo 



■ exp 



2c(t - iro) 



exp(— zciJoi ~ ikoz), 
(25) 

with fco = wo/c, for the fundamental FWM 12^]. The 
fundamental FWM is a localized, stationary free-space 
wave whose envelope propagates at luminal group ve- 
locity c, whereas the carrier oscillations back-propagate 
at the same velocity c. The eFWM is also a station- 
ary, localized wave with the same intensity distribution 
as the fundamental FWM, but propagates in a disper- 
sive medium with super- or subluminal group velocity 
l/(fcQ — a). The carrier oscillations propagate in the 
same direction at super- or subluminal phase velocity 
^o/(fco - f3)- 

Figure Efc) shows the prototype eFWM of this kind 
of wave modes, obtained from numerical integration of 
Eq. © with the approximate dispersion curve K{Q) — 
y/2koail [thick curves in Figs. Ela) and (b)] i.e., in the 
limiting case \Liu/Lp\ — 0, \Lw/Ld\ — 0), and a Gaussian 
spectrum. To pursue the validity of the model eFWM 
to describe this kind of wave modes, we have also eval- 
uated the wave mode field in some non-limiting cases 
with the same Gaussian spectrum. For \Lw/Lp\ = 1/8, 
\LwlLd\ = 1/8 [thin curves in Figs. Ela) and (b), label 1], 
the mode is nearly undistinguishable from the prototype 
eFWM, despite the dispersion curve differs significantly 



from the limiting one. Even for the relatively large ratios 
\Lw/Lp\ = 1/3, \Lw/Ld\ = 1/3 [thin curves in Figs. IHI^a) 
and (b), label 2], the calculated wave mode [see Fig. IHJd)] 
exhibits the same eFWM structure, with some incipient 
eX-wave behavior because of the actual hyperbolic form 
(not parabolic) of the dispersion curve, as explained in 
the next section. 

C. Group-velocity-dispersion-dominated case: 
Envelope X and envelope O type modes 

1. Normal group velocity dispersion: Envelope X waves 

We consider finally modes with \Ld\ <C |Lp|,|Lt„|, or 
modes of short enough duration, or propagating in a 
medium with large enough GVD. When material disper- 
sion is normal (fep > 0), the dispersion curve within the 
bandwidth approaches the X-shaped curve [see Fig. da)] 



K{n)~./hM\n\ 



(26) 



of the limiting case ji^/LpI, \Ld/Lw\ = 0. The actual dis- 
persion curve of a mode may be slightly shifted towards 
negative frequencies [as in Figs. |7|[a) and (b), labels 1 
and 2] or positive frequencies for modes with superlu- 
minal (a > 0) or subluminal {a < 0) group velocity, 
respectively. For modes with superluminal phase veloc- 
ity (/? > 0), K{^) is real everywhere [Fig. [Tfb), label 1], 
but for modes with subluminal phase velocity there is a 
narrow frequency gap about 17 = [Fig. Cfb), label 2]. A 
prototype wave mode for this case can be obtained by in- 
troducing the approximate dispersion curve of Eq. H26() 
into Eq. 0. With the two-side exponential spectrum 
f{n) = (27r/Ari)exp(-f7/Ari) we obtain 



^aA'T.T + az) ~ 3? ■ 



To 



y/kok'^r"^ + [tq + i{T + azf 



(27) 

where tq = (Ai7) ^ measures the pulse duration. Equa- 
tion H27|l is the eX wave recently described in Ref. [ijj as 
an exact, stationary and localized solution of the parax- 
ial wave equation with luminal phase and group velocities 
(a = /3 = 0) in media with normal GVD. The eX wave 
l)27|l is understood here as an approximate expression for 
modes with a,/? such that \Ld/Lp\ <C 1, |Ld/ito| *C 1- 
The spatiotemporal form of the eX wave is shown in Fig. 
He). For Ld/Lp = 1/6 (/3 > 0), Ld/L^ = 1/6 [thin 
curves in Figs. EJa) and (b), label 1], the mode retains 
an X-shaped structure [Fig. [7fd)] despite the dispersion 
curve differ significantly from the limiting one. Incipient 
PBB behavior, or radial oscillations, originates from the 
nearly horizontal dispersion curve in the central part of 
the spectrum. For Ld/Lp — —1/6 {(3 < 0), Ld/L^ = 1/6 
[thin curves in Figs. [7Ja) and (b), label 2], the X-shaped 
mode [Fig. CJd)] shows instead incipient eFWM behavior 
(light is within the cone), together with temporal oscil- 
lations arising from the frequency gap in the dispersion 
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norrmili/.ed time o 



FIG. 7; (a) Dispersion curve within the bandwidth for Ld = lO/fco, Ld/Lp 0, Ld/L^ 0, with > (thick curve), for 
Ld = lO/fco, Ld/Lp — 1/6, Ld/Lw = 1/6 (thin curve, label 1), and for Ld = lOfco, Ld/Lp = —1/6, Ld/Lw = 1/6 (thin Hue, 
label 2). (b) The same as in (a) but also outside the bandwidth of the spectrum (in arbitrary units) fiyi) — exp{—\Q\/Afl). 
(c-e) Gray-scale plots of the amplitude |?l?a,/3| of (c) the prototype eX [thick dispersion curve in (a)] with exponential spectrum 
/(fi) = exp[— (f2/A57)^], (d) of the mode with Ld/Lp — 1/6, Ld/L^ = 1/6 [thin dispersion curve 1 in (a)], and (e) of the mode 
with Ld/Lp = —1/6, Ld/Lu, — 1/6 [thin dispersion curve 2 in (a)]. Normalized coordinates are a — {t + az)AQ, p — r/ro, 
with ro = (A:ofci'An2)-i/2. 



2. Anomalous group velocity dispersion: Envelope O waves 



When \LJ « \L, 



but GVD is anomalous, 



the dispersion curve within the bandwidth can be ap- 
proached by the ellipse centered on = [Figs. Eta) 
and (b), thick curves] given by the expression 

K{n) ~ ^2ko{l3-\k'f;\ny2). (28) 

Note that the term with /3, no matter how small it is, 
must be retained to reproduce the real-valued part of the 
dispersion curve. The group velocity of the mode can be 
shghtly subluminal {a < 0) or superluminal {a > 0), as 
in Fig. Eta) and (b) (thin curves), but the phase veloc- 
ity of these modes is always superluminal {(3 > 0). An 
approximate analytical expression for this type of modes 
can be obtained by introducing the approximate disper- 
sion curve of Eq. (|28|l into Eq. (|HJ|. Under condition 



l^rfl ^ l^pl, the frequency gap fig ~ ^2/3/|fcQ | is much 
smaller than Af2, so that the amplitude spectrum /(f2) 
can be assumed to take a constant value in the integra- 
tion domain of integral in Eq. (jS]) , which then yields the 
expression 



1 



X sm 



^2p/\k'^\^ko\k'^\r^ + {T + az) 



, (29) 



of the same form as the 0-type impulse response mode 
in media with anomalous dispersion. Figure EJc) shows 
its spatiotemporal form. For comparison, the wave mode 
with Ld/Lp = -1/6, Ld/Ly, = -1/8 [Fig. Ela), thin 
curve] and the two-sided exponential spectrum [Fig. |S| 
(b)] was calculated from Eq. ||SJ), and its 0-shaped spa- 
tiotemporal form is depicted in Fig. Eld). 



IV. NONPARAXIAL DESCRIPTIONS OF WAVE terms of the characteristic lengths remains essentially un- 
MODES 

The purpose of this section is to show that the preced- 
ing classification of wave modes in dispersive media in 
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(c) 



-S 5 

lUHinaliziMl lime 



nonnalizcd lime 6 




FIG. 8: (a) Dispersion curve within the bandwidth for |Ld|/|Lp| —> 0, — > 0, with Ld < (thick hne), and for 
Ld/Lp = —1/6, Ld/Lw = —1/8 (thin curve), (b) The same as in (a) but also outside the bandwidth of the spectrum (in 
arbitrary units) f(yi) — exp(— jfij/Afi). (c) Gray-scale plot of the amplitude I'&a./gl of the eO wave of Eq. 1291 . (d) and of 
the mode with Ld/Lp — —1/6, Ld/L^ = —1/8 and the exponential spectrum of (b), numerically calculated from Eq. (|HJ. 
Normahzed coordinates are a = [r -\- az) y'2/3/|fc(/| , p = \/2kol3r. 



altered when performed from the more exact nonparax- 
ial approach, if condition jSJ of quasi-monochromaticity, 
and and Q of quasi-luminal group and phase veloc- 
ities are satisfied. 

We consider now the polychromatic Bessel beam, 

E{r,z,t) = 77- / (Ljf{uj~uJo) 
JK real 

X Jo{Kr) exp{ikzz) exp{—iu;t), (30) 

where K and must be related hy K = yjk'^iuj) — for 
each monochromatic Bessel beam component to satisfy 
the Helmholtz equation AE + k'^{u))E = 0. Stationarity 
of the intensity in some moving reference frame requires 
the axial propagation constant kz to be a linear function 
of frequency Q , a condition that is suitably expressed as 

k,{n) = {ko- P) + {k'^~a)n. (31) 

Equation (|30() can be then rewritten in the form 
E{r, z, t) = 3>a,/3(?'7 T + Oiz) exp{—iPz) exp(— iwot + ik^z), 
where the reduced envelope is given by the same expres- 
sion as in the paraxial case, namely, 

$„,^(r,r + az) = / dnf{n) 
2tt JK{n) real 

X Jo[ii:(r2)r]exp[-ir2(T + az)],(32) 

but with a transversal dispersion relation K{^) = 
\J k^{il) — A:^(f7) given now by 

K{n) = [{2koP- (3'^) + 2{koa + k'o(3-af3)Q 

+ {kok'^ + 2k'oa-a^)n^Y^^ (33) 



up to second order in dispersion [k(fl) = ka + k'gQ + 
k'^n^/2]. 

In the case of propagation in free-space (/cq — ^o/c, 
/cq — 1/c, /cq = 0, with c the speed of light in vacuum), 
Eqs. 123 and (ESJ yield Eq. (7) of Ref. for general 
free-space FWMs, if the identifications a = (1 — 7)/c and 
P — ivoa + 2jl3s are made (7 and /3s being the the free 
parameters defined in Ref. 23). In particular, the case 
with a = yields the original Brittigham's FWM |2 l l2 ^ . 
and the case with (3 = ojQa yields the Bessel-X pulse of 
cone angle 9 = (2ca)^/^, or X wave with narrow spectral 
amplitude centered at an optical frequency, introduced 
by Saari in Ref. 16, and demonstrated in Ref. IT^ 

In a dispersive media, and under conditions I^J and 
^ of quasi-luminality, we can neg lect in Eq. ^ the 
terms (3"^, afi and in comparison with 2fco/3, k^a and 
2kQa, respectively, to obtain the approximate expression 

K{n) ~ ^2{ko + k'„n) {13 + an) + kQk'^9? (34) 

for the nonparaxial dispersion relation of quasi-luminal 
modes. The first conclusion is then that the paraxial 
dispersion curve [Eq. Q] may significantly differ from 
the nonparaxial one [Eq. H34(l ]. even if conditions Q 
and ^ are satisfied. In fact, it is not difficult to find 
set of parameters for which the nonparaxial dispersion 
curve is, for instance, a vertical hyperbola, whereas the 
paraxial dispersion curve is an horizontal hyperbola [see 
Fig. Ha)]. 

From a physical point of view, however, it is only the 
portion of the dispersion curve within the mode band- 
width that is of relevance for the spatiotemporal mode 
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the (exactly) horizontal straight line 



FIG. 9: (a) Paraxial and nonparaxial transversal dispersion 
curves of the modes of carrier frequency ujq = 4 fs~ 1 with a — 
300 mm~^fs and /3 = 400 mm~^ in fused silica (fco — 19530 
mm-\ k'o = 4988 mm-^fs and k'o = 77 mm-^fs^). (b) The 
same as in (a) but only within the bandwidth of the shortest 
(widest spectrum), single-cycle wave mode (AO = ujo/2tt). 



structure, and, as our second conclusion, this portion is 
approximately the same in the paraxial and nonparax- 
ial approaches if the additional condition ^ of quasi- 
monochromaticity is also satisfied: Writing, for trans- 
parent dispersive materials, /co/Zcq ~ ^Oi we obtain 



K{n) « ^j2ko (1 + n/coo) (/3 + an) + kokl^n^, (35) 
or, in terms of the mode characteristic lengths. 



'2fco 



An 



1+ — n,, {Lp 



UjQ 



Lw nr. 



(36) 

Since (Ari/ti;o)|f^n| <C 1, the nonparaxial dispersion 
curve within the bandwidth can be approached by the 
paraxial one, that is, by Eq. H19|l . as illustrated in 
FigEIb) for the extreme case (widest possible bandwidth) 
of a single-cycle mode (Afi/wo = 1/2-11). In particular, we 
can afhrm that the description performed in Section IlIII 
of quasi-monochromatic, quasi-luminal modes in terms 
of their characteristic lengths is independent of the ap- 
proach used. 

To illustrate the relationship between the paraxial and 
nonparaxial approaches, and the type of results we can 
expect from the paraxial one, we consider wave modes of 
any bandwidth AS! propagating in normally dispersive 
media [k'^ > 0) with 



- k' 



jL/2 
ftp 



{k'o - + koK 
~7W^k^ 



(37) 

kXfik'o^ (38) 



[see Figs. llOf al and (b) for propagation in fused silica] , so 
that the nonparaxial dispersion curve is, from Eq. (|34|l . 



K{n) = K = 



i,3 I." 



k'o kok^ 



(39) 



and the corresponding nonparaxial wave modes are the 
dispersion- free, diffraction-free PBBs $c(,/3(''', t + az) = 
f{T + az)Jo{Kr) studied in Ref. 10]. The approximate 
equalities in Eqs. (EHl and hold for weakly 

dispersive materials such that fcg <C k'^ /k^, in which case 
a and /? satisfy conditions Q and ((SJ of quasi-luminality 
for the group and phase velocities. As seen in Figs. llOf a') 
and (b), this is the case of fused silica at any visible 
carrier frequency. 

For these PBBs, it is easy to see that the paraxial 
and nonparaxial descriptions become undistinguishable, 
in spite of the apparent drawback that PBBs are no 
longer exact solutions of the paraxial wave equation in 
dispersive media [when fcp ^ 0, the paraxial dispersion 
curve ^ is never an horizontal straight line]. In fact, 
when fcg <C k'^/ko, the relationship <C ^ l^dl 
is satisfied for any mode bandwidth down to the single- 
cycle limit [see Fig. llOf c) for the case of fused silica]. 
Accordingly, these modes are of PBB type, that is, the 
paraxial dispersion curve within the bandwidth can be 
approached by an horizontal straight line [see Fig. Iior dl 
for luq = 2 fs~^ in fused silica]. Finally, the paraxial pro- 
totype PBB for these modes is given, from Eq. (|21() . by 
<l>aAr, r+az) = f{T+az)Jo{Kr), with K = ^/k^^Jk^, 
that is, by the same expression as in the nonparaxial ap- 
proach. 



V. CONCLUSIONS 

Summarizing, we have described and classified the 
pulsed versions of Bessel beams with the property of be- 
ing localized and remaining stationary (diffraction-free 
and dispersion-free) during propagation in a dispersive 
material with slightly super- or subluminal phase and 
group velocities. As for the wave mode description, 
we have found the analysis of the transversal dispersion 
curve K{n) to be an useful tool to understand the spa- 
tiotemporal mode structure. Wave modes have been clas- 
sified into three broad categories: PBB-like, eFWM-like, 
and eX-like (eO-like) modes, depending on the relative 
strength of their phase and group velocity mismatch with 
respect to a plane pulse, and defeated GVD, as measured 
by the mode phase- mismatch length Lp, group- mismatch 
length and the dispersion length Ld- 

We have verified that the paraxial description leads 
to the same description and classification as would be 
obtained from the more accurate nonparaxial approach 
when the conditions of narrow bandwidth (j^J and of 
quasi-luminality H4I5() are satisfied. All previously re- 
ported optical Bessel beams, X waves, Bessel-X waves, 
or focus wave modes generated by linear or nonlinear 
means satisfy indeed these requirements. 
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paraxial dispersion curves within the bandwidth, given, re- 
spectively, by Eqs. and 13311 . 
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